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Constituent quark scaling violation due to baryon number transport 

J.C. Dunlop,' M.A. Lisa,^ and P. Sorensen' 

Brookhaven National Laboratory, Upton New York 11973, USA 
^Department of Physics, Ohio State University, Columbus, Ohio 43210, USA 

In ultra-relativistic heavy ion collisions at ,/snn ~ 200 GeV, the azimuthal emission anisotropy of hadrons 
with low and intermediate transverse momentum (pj < 4 GeV/c) displays an intriguing scaling. In particular, 
the baryon (meson) emission patterns are consistent with a scenario in which a bulk medium of flowing quarks 
coalesces into three-quark (two-quark) "bags." While a full understanding of this number of constituent quark 
(NCQ) scaling remains elusive, it is suggestive of a thermalized bulk system characterized by colored dynamical 
degrees of freedom- a quark-gluon plasma (QGP). In this scenario, one expects the scaling to break down as the 
central energy density is reduced below the QGP formation threshold; for this reason, NCQ-scaling violation 
searches are of interest in the energy scan program at the Relativistic Heavy Ion Collider (RHIC). However, as 
^/sNN is reduced, it is not only the initial energy density that changes; there is also an increase in the net baryon 
number at midrapidity, as stopping transports entrance-channel partons to midrapidity. This phenomenon can 
result in violations of simple NCQ scaling. Still in the context of the quark coalescence model, we describe a 
specific pattern for the break-down of the scaling that includes different flow strengths for particles and their 
anti-partners. Related complications in the search for recently suggested exotic phenomena are also discussed. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq 

Keywords: quark-gluon plasma, constituent quark scaling, stopping, heavy ions 



I. INTRODUCTION 



Quark coalescence in the highest-energy heavy ion 
collisions 



Hadronization- the process through which a state charac- 
terized by colored dynamical partons is resolved into a state 
of colorless hadrons- is central to the theory of the Strong 
interaction, but remains only incompletely understood. An 
important aspect of this process has become clear in studies 
of the forward region in hadron-hadron collisions |1| and in 
high-multiplicity collisions of ultra-relativistic heavy ion col- 
lisions |l2][3l: the hadronization of a parton is strongly af- 
fected by the presence of other partons close in phase space. 
Whereas the vacuum hadronization of a single paiton liberated 
in a high momentum-transfer (Q^) interaction is described in 
terms of string-breaking scenarios or parameterized in frag- 
mentation functions, there is mounting evidence that, in a 
dense phase-space scenario, colored partons essentially "coa- 
lesce" into colorless bound states, much like the formation of 
light nuclei (e.g. deuteron or triton) from free nucleons emit- 
ted from a hot zone |4-6|. Models based on this remarkably 
simple mechanism, not understood at a fundamental level, 
have enjoyed considerable success at describing the "leading 
hadron effect" Qj as well as the multiplicity dependence of 
yields, spectra and momentum anisotropics from heavy ion 
collisions at the highest energies at the Relativistic Heavy Ion 
Collider (RHIC) 



non-valence quarks. Indeed, the valence quarks^- three for 
baryons and two for mesons- were originally invented to ex- 
plain the flavored quantum numbers- isospin, strangeness, 
etc- of the hadrons. Only later was the connection established 
between the valence quarks and the high momentum-fraction 
(x) fermionic partons. Nevertheless, the two or three valence 
quarks represent the lowest Fock states of the partonic wave- 
function of a hadron, and these appear the relevant degrees 
of freedom; it is argued that inclusion of higher-order Fock 
states does not significantly modify the description of the co- 
alescence process and the related phenomenology 1 13|. 

The data at RHIC is consistent with a partially thermal- 
ized system of deconfined quarks undergoing collective ex- 
pansion with an azimuthal anisotropy in momentum space 
proportional to the initial spatial anisotropy 13] [l4l [TT| . As 
the system cools, pairs and triplets of neighboring quarks co- 
alesce to become the valence quarks (or "constituent quarks") 
of mesons and baryons (where the gluons presumably con- 
tribute to dressing the valence quarks |16 17 1). The origi- 
nal flow pattern of the deconfined quarks leaves a simple and 
characteristic fingerprint on the momentum distribution of the 
observed hadrons, since a hadron's momentum is simply the 
(vector) sum of the momenta of its valence quarks 
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where n — 2 (3) for mesons (baryons). In the simplest instan- 
taneous 2 ^ 1 or 3 — > 1 coalescence process, only three of the 
four momentum components are conserved; either energy or 
momentum conservation is violated lfT2l . The most important 



Even more remarkably, the objects that coalesce appear 
to be valence quarks. At first, this is surprising, since the 
dynamical quantities in QCD (i.e. the ones that caiTy mo- 
mentum) are partons which are overwhelmingly gluons and 



' Anti-quarks are treated on equal footing with quarks in coalescence mod- 
els. Dynamically, we will treat antiquarks as simply another variety of 
quark. 



features of quark coalescence that we discuss in this paper are 
not substantially altered in a more complete treatment of the 
phenomenon, accounting for energy and entropy conservation 
effects 1 12;'T8', 191. 

In particular, the quarks' radial flow- the enhancement to- 
wards higher transverse momentum (pr) due to pressure- 
driven expansion of the bulk source- is reflected most strongly 
in the three-quark baryons. Thus, coalescence provides a nat- 
ural explanation for the "anomolous baryon enhancement" 
at intermediate pr observed in the highest-energy heavy ion 
collisions at RHIC ^. Similarly, the nuclear modification 
factor- the scaled ratio Raa (pt) of transverse momentum dis- 
tributions from heavy ion and p + p collisions- shows a clear 
separation into mesons and baryons ||3] \T5\ . 

In heavy ion collisions, the azimuthal anisotropy of the mo- 
mentum distribution is characterized by Fourier components. 



— -oc 1+2 > v„(;?7-)cosn(|), 



(2) 



where (|) is measured relative to the direction of the impact pa- 
rameter [c.f.|20l for a full discussion]. Of particular interest is 
the elliptic flow parameter vj (pt), which is strongly sensitive 
to the equation of state of QCD matter (e.g. speed of sound) 
as well as tranport coefficients like viscosity. For small values 
of V2 and narrow hadronic wavefunctions, the elliptic flow pa- 
rameters of a bulk system of quarks (a and b) and the mesons 
into which they coalesce, are related by lfT2ll 



vf (pr) = V2 [xaPr) + v'2 {xbPT) 



(3) 



for fixed momentum fractions Xa and X), {xa +Xfc = 1), with an 
analogous equation for baryons. 

In the event that the constituent quarks (a and b) have the 
same elliptic flow before hadronization. Equation [3] leads to 
the number-of-constituent-quark (NCQ) scaling pattern ob- 
served at RHIC, 



v^l 



p't 



{Pr/n 



(4) 



where n — 2 (3) for mesons (baryons). In this scenario, all 
mesons should follow one common V2{pt) curve, and all 
baryons another. The two should be related via 



vf(/^r/3) vfipT/2) 



^Upt)). 



(5) 



B. Violation of NCQ scaling 

The observed satisfaction of Eq. l5] in AuH-Au collisions at 
RHIC was one of the most compelling indications that de- 
confined partonic degrees of freedom were playing a dynamic 
role in the bulk medium of the early phase. Consequently, the 
breakdown of this scaling as the beam energy is reduced, is 
an important signal in the RHIC energy scan program fT\\. 
It may indicate that the initial energy density of the system is 
below the threshold for QGP creation, pinpointing the phase 
transition between confined and deconfined QCD matter 



Additionally, it has been recently proposed f22l that a de- 
confined bulk system at finite baryon density may acquire an 
electric quadrupole moment due to chiral magnetic waves in 
the plasma. This interesting phenomenon would break the 
degeneracy between V2 for positive and negative pions, also 
clearly breaking NCQ scaling. 

It is important to consider other less exotic mechanisms 
that may also cause violation of the scaling represented by 
Eq. Is] In particular, we recall that this scaling should hold 
if all quarks (and antiquarks) have the same underlying flow 
distribution (vj in equation Db. This would be a natural con- 
seqence of thermalization. The scenario where all constituent 
quarks have the same Vj we call NCQi . 

In this paper, we examine whether the breakdown of Eq. |5] 
would necessarily signal that the hadrons are not resulting 
from the coalescence of flowing constituent quarks. In par- 
ticular, we discuss a minimal extension of the NCQi model 
that retains constituent quarks as the dynamical degrees of 
freedom and coalescence as the hadronization mechanism. 
However, the assumption that all quarks have the same vi^ is 
discarded due to the well-recognized phenomenon of baryon 
stopping, which is increasingly important at lower energies. 
In particular, we recall that the transport of baryon number 
from the entrance channel to midrapidity ("stopping") is in- 
creasingly important at lower energies. Since we continue to 
work in the dynamical constituent quark paradigm, this means 
that u and d quarks transported from y = jbeam to y = surely 
suffer multiple collisions with each other. Meanwhile, at the 
lower energies in question, the quark-antiquark pairs created 
from the vacuum may experience relatively fewer In a pic- 
ture where Vj is developed through collisions, it is not unrea- 
sonable to expect that quarks transported from y = jbeam to 
y = Q will develop a larger V2- 

Such considerations will lead to a specific pattern for the 
breakdown of Eq.jS] Without invoking exotic phenomena, this 
simple scenario also implies that the degeneracy of V2 for par- 
ticles and their anti-partners will be broken in a specific way. 
Not only its prediction for V2 [7i+] — V2 [ti^], but also that for 
V2 [K^] — v'2 \K^], can be compared to data and predictions 
from more complicated models. 

In the following section, we briefly review the energy de- 
pendence of stopping in heavy ion collisions- transport of 
baryon number from the high-rapidity region of the entrance 
channel to midrapidity in the exit channel. We also use 
measured particle yields to estimate the fraction of u and d 
quarks at midrapidity that would arise from baryon stopping 
at two collision energies. In section III we consider quan- 



titatively a two component model for quark number scaling 
NCQ2, in which the phenomenon of stopping leads to at least 
two samples of quarks with different Vj values which then co- 
alesce into hadrons. For tractability, we idealize this non- 
thermal distribution in a two-component formalism: trans- 
ported quarks follow one flow profile and produced quarks 



another. We briefly summarize in section IV 



n. TRANSPORT OF ENTRANCE-CHANNEL QUARKS TO 
MIDRAPIDITY 



In this section, we briefly review the phenomenon of baryon 
stopping in heavy ion colHsions. We then use hadron yields 
measured by the NA49/SPS coflaboration to estimate the frac- 
tion of u and d quarks at midrapidity that would arise from 
stopping. 



These fractions will be used in section III 



as mput 

to a simple model to predict the breakdown of NCQ scaling, 
Eq.|5] 



A. Stopping in heavy Ion collisions 

In the RHIC energy scan program, the beam energy is var- 
ied to modify the initial conditions of the hot QCD system 
created. In addition to changing the energy density of the ini- 
tial state, it is well-accepted that, due to baryon stopping, the 
baryo-chemical potential /jb of the system is larger at lower 

Baryon stopping, the transport of baryon number from its 
initial location at beam rapidity towards the initially baryon- 
free region at mid-rapidity, is most directly measured via the 
rapidity distribution of net protons (the number of protons mi- 
nus antiprotons). At ,/snn ~ 5 GeV, the rapidity distribution 
is peaked at midrapidity. As the collision energy is increased, 
the distribution peaks at increasingly forward rapidity. This 
behavior has been parametrized as an average rapidity loss, 
which increases from approximately 1 unit at y^imv ~ 5 GeV 
to 1.7 unit at ^/snn ~ 17 GeV, with a smaller rise towards 2 
units by the highest RHIC ^/sj^ of 200 GeV||23]. This rise 
in rapidity loss is less rapid than the rise in the beam rapidity 
with increasing y^s^N, leading to a decreasing population of 
net baryon number at mid-rapidity. 

More detailed statistical model calculations, based on mea- 
surements of the yields of a range of particles, agree with 
this general conclusion. Within these models, with increasing 
y/sjsiN the net baryon density first rises, achieves a maximum 
at y^SNN w 8 GeV, and then falls for higher ^/snn ll24l . The 
system transitions from one with entropy density dominated 
by baryons at low energy to one dominated by mesons at high 
energy, with equal fractions at y^s^N ~ 8 GeV 1251 . The frac- 
tional importance of transported quarks to the system's evo- 
lution grows rapidly with decreasing y^s^N in the lower end 
of the region probed by the RHIC Beam Energy scan, from 
y^s]^ w 15 GeV to 7.7 GeV. 



Estimates for the stopping contribution to Ught quark 
yields at midrapidity 



In section III we develop a simple model in which trans- 



ported u and d quarks have stronger flow than do produced 
quarks (including produced u and d quarks). An important in- 
gredient to this model is the fraction of u (d) quarks present 
at midrapidity that arise from baryon number transport. In 
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TZ+ 


72.9 


72.9 








72.9 




%^ 


84.8 




84.8 




84.8 






n"(*) 


78.85 


39.43 


39.43 




39.43 


39.43 
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1.17 






1.17 






1.17 


K+ 


16.4 


164 










164 


K- 


5.58 






5.58 


5.58 






K0(*) 


10.99 




10.99 








13.84 


K0(*) 


10.99 






10.99 




10.99 
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46.1 


92.2 


46.1 










n(*) 


70.84 


70.84 


141.68 
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13.4 


134 


13.4 


13.4 








E^ 


0.93 




0.93 


1.86 








P 


0.06 








0.12 


0.06 




n(*) 


0.04 








0.04 


0.08 




A 


0.1 








0.1 


0.1 


0.1 


Sum 




305.17 


337.33 


33 


130.07 


123.56 


28.66 



TABLE I: Left two columns: midrapidity yields of com- 
mon particles from central Pb+Pb collisions measured by the 
NA49/SPS Collaboration |26-33| at ^/S^ = 6.41 GeV. Starred 
hadrons are not measured, but estimated from other hadrons. 
In particular, dN[lt^^]/dy = 0.5{dN[ll+] /dy + dN[K-] /dy), 
dN [k"] /dy = dN [K°] /dy = 0.5 [dN [K+] /dy + dN [K ] /dy), 
dN[n] = 1.54dA'[p], and dN[n] = \.5A-'^dN\p]. The factor 1.54 is 
the neutron-to-proton ratio of Pb. Right six columns: midrapidity 
yield of constituent quarks for each hadron. 



particular, we want the fraction 



X„T — 






(6) 



where N^,t is the number of u quarks from the incoming heavy 
ions transported to midrapidity, and A^„/> is the number of u 
quarks produced from u — il pair production at midrapidity. 
The fraction X^t is defined similarly. 

To estimate X^j and X^t, we use measured midrapidity 
yields of common particles from central Pb-nPb collisions by 
the NA49/SPS Collaboration at y^s/^ = 6.41 and 8.86 GeV. 
Tables |l] and |ll] list the measured yields of hadrons and their 
constituent quarks. 

Since produced quarks and antiquarks are formed in pairs, 
the number of transported u quarks is given by the imbal- 
ance between the total number of u and m quarks; X^r — 
{Nu—Nu) /Nu- Table III lists the fractions for the two ener- 
gies. 

Several comments are in order. Firstly, at these energies, 
roughly half of the light constituent quarks at midrapidity 
orginate from the colliding nuclei; clearly stopping cannot be 
ignored. Secondly, the fraction of d quarks transported from 
the y = jbeam is greater than the fraction of u quarks, simply as 
a consequence of the isospin of the entrance channel. Thirdly, 
the imbalance in s and s quarks in tables [I] and [ll] reminds us 
that our estimates are just that. 
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13.84 
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13.84 
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13.84 
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82.6 


41.3 
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63.5 


63.5 


127 
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1.02 


3.06 


3.06 


3.06 








A 


14.6 


14.6 


14.6 


14.6 








E^ 


1.15 




1.15 


2.3 
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0.32 








0.64 


0.32 




n(*) 


0.21 








0.21 


0.42 




A 


0.33 








0.33 


0.33 


0.33 


S+ 


0.07 










0.07 


0.14 


Sum 




331.14 


357.73 


39.48 


165.54 


162.26 


35.57 



TABLE II: The same as for tablell] but for 



'SNN ' 



8.86 GeV colli- 



\/SNN 


X„7- 


XjT 


6.41 GeV 
8.86 GeV 


0.57 
0.50 


0.63 
0.55 



TABLE III: Based on the data in tables|l]and|ll] the fraction of u and 
d quarks at midrapidity, that originate from stopping of quarks in the 
colliding Pb nuclei. See text for details. 



III. VIOLATIONS OF SIMPLE NCQ SCALING IN A 
TWO-COMPONENT SCENARIO 

As discussed in section IT] the violation of NCQ scaling, 
eq.|5]and the breakdown of the degeneracy of V2 between par- 
ticles and antiparticles, implies a lack of kinetic thermaliza- 
tion of the dynamic system. This may arise if the energy of the 
collision falls below the threshold to produce a flowing system 
of deconfined partons, so that particle-specific hadronic cross- 
sections determine each hadron's flow strength. Alternatively, 
it may arise from more interesting phenomena such as chiral 
magnetic waves. Our approach in this paper is a minimalist 
one, asking whether such effects may be expected without in- 
voking exotic phenomena or abandoning a scenario of flowing 
quarks which coalesce into flowing hadrons. 

The absence of complete kinetic thermalization at low en- 
ergies (y^ivviv ^ 30 GeV) has long been recognized, based on 
the severe failure of hydrodynamics to reproduce hadronic 
V2 |34|. In short, at these energies, the constituents of the 
system do not rescatter sufficiently to achieve thermalization. 
As we have discussed, at these energies baryon transport from 
the entrance channel plays a huge role. We wish to test the ro- 
bustness of the constituent quark paradigm, so the transported 
baryon number is represented by transported u and d quarks. 



The very fact that they have been transported over a significant 
rapidity range attests to the likelihood that these quarks, in any 
event, have suffered many scatterings. We make the plausible 
postulate that transported quarks experience more scatterings 
than produced ones at these energies, hence approaching the 
thermal limit more closely and developing a larger V2- 

Clearly, the resulting non-thermal quark momentum distri- 
bution reflects a continuum of quarks rescattering more or less 
before coalescence. In order to render the problem tractable, 
we model the situation in a simple limit of two populations: 
transported quarks and produced quarks, with the former pop- 
ulation characterized by a stronger flow than the latter We 
emphasize that this is a simplification in order to make a point: 
we do not propose that there are really two distinct thermal- 
ized fluids created in a heavy ion collision. 

Hence, we have two populations of constituent quarks with 

FT _ 

distinct flow fields, Vj and V2 , for produced (m^, m^, t/^, c/^, 
5^, s'') and transported (m^, d^) quarks, respectively. In this 
simplest two-component model, a hadron's elliptic flow pa- 
rameter is given by 



4 (Pt) = t kf ^'2' (PT/n) + (1 -X^j) 4 {PT/n) 



, (7) 



where X t is the fraction of quark species o, that originates 
from baryon stopping, as discussed in section O As per the 
discussion in that section, reasonable estimates are X„ = 0.50 
and Xd = 0.55. NaturaUy, X,^r = Xjt = X,r = X^t = 0. 

Figure [T] shows a example of the resulting V2 from our sim- 
ple NCQ2 scenario. For the purpose of illustration, for these 
calculations, we had to assume some functional form for the 
quark elliptic flow. We chose the same functional form for 
both produced and transported quarks: 



V2 [pt) = Mtanh {pj / (0.5 GeV/c)) . 



(8) 



For the example in figure [T] M = 0.07 for transported quarks 
and M ~ 0.05 for produced quarks. The choice of this partic- 
ular functional form is rather arbitrary and does not affect the 
points we make below. 

Clearly, simple NCQ scaling (equation [5]) is violated, with 
an unavoidable species-dependent signature. In particular, 
one finds: 

V2 [kT — du\ > V2 \ji^ = ud\ 

V2 [K^ — us\ > V2 [K^ = us\ 

V2[p ~ uud] > \'2 \p = uucl\ (9) 

V2 [A — uds] > V2 [a = uds\ 

V2[p — uud] > V2 [A = uds] 

(v2 [p — uud] ~V2[p — uud] ) > (v2 [A = uds] — \'2 [A ~ udf^ ) , 

It is interesting that the ordering of V2 for positive and nega- 
tive pions is the same as that predicted due to chiral magnetic 
wave effects |22|. We also find that the charge-ordering for 
kaons (v2 [K+] > V2 [K^]) is opposite to that for pions. The 
chiral magnetic wave effect would generate the same charge- 
ordering for pions as for kaons, thus providing a testable dis- 
tinction between the chiral magnetic model and our stopping- 
based model. However, hadronic effects (e.g. the smaller 
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FIG. 1: (Color online) NCQ2, the simplest generalization of the 
NCQi model, in which the transported up and down quarks have a 
40% stronger intrinsic V2 than do the produced quarks, which them- 
selves all have the same V2 . The fraction of u and d quarks that are 
transported is 50% and 55%, respectively. See text for details. 
Panel (a): Intrinsic V2 of all quarks. 

Panel (b): Hadron V2 based on coalescence of quarks shown to the 
left. 

Panel (c): NCQ-scaled hadron flow, V2 [pr/n) jn. 
Panel (d): Deviation from simple NCQ scaling- the curves from the 
lower left panel, divided by their average. 

In panels (b)-(d), K+ data points lie beneath those for n+, K^ data 
points lie beneath those of (|), and the data points for all anti-baryons 
are coincident. 



cross-section for K ) may complicate the interpretation of 
such a test |3.5j. 

Comparing the anisotropics of particles and their anti- 
partners, as listed in equation |9j is straightforward and rela- 
tively unambiguous. The details of cross-species comparisons 
can depend more on the particular functional forms used for 
the quark flow profiles, the weighting factors Xj and whether 



the latter depend on pj. For the simple case we have con 

2 



sidered, ^^ 



> 



; the entire species dependences can be 
seen in panels (d) of figures [T] and |2] 

Our primary points have been made already in this simple 
model, but we mention an additional complication. If there is 
insufficient rescattering to fully thermalize the light produced 
quarks, then the heavier strange produced quarks are likely to 

be even less thermalized. In this case, v^ < Vj ''^ ; similar 
considerations have been discussed by Lin and Ko f361. Fig- 
ure [2]shows the situation when the functional form of equa- 
tion [8] describes the flow of all quarks, as before, but now 
M=0.045 for the strange quarks. In this case, the degeneracies 



FIG. 2: (Color online) The same as in figure [T] but the produced 
strange quarks have 10% less intrinsic V2 than do the produced light 
quarks. See text for details. 



(e.g. V2 [K^] = V2 [(])]) seen in figure [Tland listed in its caption 
are broken; all hadrons have distinct elliptic flow curves. 

Additional reasonable complications can be considered. 
Clearly, the functional forms used for the quark flow can be 
varied from the simple form (eq. [8]l used here. Furthermore, 
one may reasonably argue that the fraction of light quarks aris- 
ing from transport {X t) should depend on pj', we have treated 
it as a constant for simplicity. Exploring such considerations 
amounts to tuning the model. We leave such explorations for 
later comparison and fitting when data become available. 



IV. DISCUSSION AND SUMMARY 

The success of NCQ scaling of elliptic flow at 
^snn=2QQ GeV has been one of the most striking observa- 
tions at RHIC, strongly suggesting the creation of a flowing, 
thermalized bulk system of quarks that coalesce into hadrons. 
Hence, observing the violation of this scaling as ^Jsnn is de- 
creased could be of crucial importance, both for validating 
the simple dynamical constituent quark model, and for pin- 
pointing the conditions required to undergo the deconfinement 
phase transition. Furthermore, recent theoretical predictions 
suggest that a chiral magnetic wave effect may reveal itself by 
inducing a different flow for positive and negative pions |22|. 
The observation of NCQ scaling violations would thus be po- 
tentially far-reaching. 

It is important, therefore, to explore less exotic reasons for 
any scaling violations. We have discussed one simple scenario 
here, which requires neither a fundamental difference in the 



phase of QCD matter in the measured energy range, nor a new 
exotic effect. 

The model predicts an unavoidable species-dependent pat- 
tern for the breakdown of NCQ scaling and depends on only 
two assumptions. Firstly, it assumes that, just as at top RHIC 
energies, the system can be described in terms of constituent 
quarks that coalesce into hadrons as the system cools. Sec- 
ondly, it assumes that quarks transported from beam rapidity 
to midrapidity suffer more violent scatterings than do quarks 
produced at midrapidity at low y^s^N- We simplified the situ- 
ation by treating the system as two distinct quark populations, 
but our main points do not depend on this simplification. 

(Baryon transport from the entrance channel is another im- 
portant ingredient of the model, but its relevance is far from an 
assumption; the phenomenon of stopping is well known and 
the isospin effect (X^t > X^t ) is obvious and based on data, 
as discussed in section [H]) 

The second of our two assumptions seems at least very 
plausible. It is clear that at low energies, the system does 
not have sufficient density or energy to fully thermalize- the 
dynamical constituents do not scatter enough. Unlike the pro- 
duced particles born at midrapidity, however, the transported 
quarks had to undergo several collisions just to reach midra- 
pidity, after which they could rescatter further 



It is the first assumption- that even at low energies 
where scaling violations might be found, the system is well- 
described by a flowing system of constituent quarks- that 
seems most questionable. Nevertheless, our task has been 
to explore the implications of its validity even at low ^snn- 
We have found an unambiguous species dependendence of 
V2 listed in equations [9] Quantitative details depend on tun- 
ing which we do not consider in this first study. Detailed 
comparisons with experimental data should be performed, but 
we have shown that violation of NCQ scaling or particle- 
antiparticle v'2 degeneracy themselves is insufficient to claim 
either the crossing of the deconfinement threshold or exotic 
phenomena. 
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